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We present a stability theory for kink propagation in chains of coupled oscillators and a new 
algorithm for the numerical study of kink dynamics. The numerical solutions are computed using 
an equivalent integral equation instead of a system of differential equations. This avoids uncertainty 
about the impact of artificial boundary conditions and discretization in time. Stability results also 
follow from the integral version. Stable kinks have a monotone leading edge and move with a velocity 
larger than a critical value which depends on the damping strength. 

PACS numbers: 05.45.-a; 83.60.Uv; 45.05.+X 



I. INTRODUCTION 

The dynamics of waves in chains of coupled oscilla- 
tors is the key to understanding the motion of defects in 
many physical and biological problems: motion of dislo- 
cations jl 3l or cracks [j in crystalline materials, atoms 

motion of electric 



adsorbed on a periodic substrate 
field domains and domain walls in semiconductor super- 



lattices |j| , pulse propagation through myelinated nerves 
[^| or cardiac cells [?J... A peculiar feature of these spa- 
tially discrete systems is that wave fronts and pulses get 
pinned for entire intervals of a control parameter such 
as an external force. Typically, wave fronts do not move 
unless the external force surpasses a control value. Such 
is the case with the static and dynamic Peierls stresses 
in dislocation dynamics |^| or the dynamic and static 
friction coefficients j^| in continuum mechanics. Pinning 
and motion of wave fronts also explain the relocation of 



static electric field domains and the self-oscillations of 
the current in semiconductor superlattices |5|. 

Wave front motion in systems of nonlinear oscillators 
modeling these phenomena are easier to analyze in the 
overdamped case, and less so if inertia is important. In 
the presence of inertia, the wrong choice of boundary 
conditions or the numerical method may suppress im- 
portant solutions of the original system or yield spurious 
oscillations. Thus two problems that are important in all 
spatially discrete systems acquire even more importance: 
how do we find wave fronts? and, what are their stability 
properties? 

We have solved the first problem in a recent work 
by choosing a damped system of oscillators with a piece- 
wise linear source term, see also 3, 11 , ^ . Our results 
show explicitly the existence of kinks with oscillatory pro- 
files for systems with little or no damping. In the latter 
case, these wave fronts have at least one tail with non 



decaying oscillations that extend to infinity. Depending 
on the control parameter, branches of oscillatory wave 
fronts may exist, coexisting for entire intervals of the ex- 
ternal force and even coexisting with pinned wave front 
solutions. These facts, long-lived oscillatory profiles and 
coexistence of wave front branches, highlight the impor- 
tance of ascertaining the stability properties of these so- 
lutions. This is not easy and not many results are known. 
To be precise, let us consider oscillator chain: 

mu'^+au' n = K{u n+ i - 2u n + u n -i) — V'(u n )+W. (1) 

We nondimensionalize the model by choosing the time 



scale y , where a is the interatomic distance, m is 
the mass and v m the strength of the on-site potential. 
For a piecewise parabolic potential, the nondimensional 
equation is: 

u n + l u 'n = D{u n+ l - 2u n + Un-l) - ff(«n) + F, (2) 

1 



s s < 



s-1 s> 



(3) 



Here 7 and D are the ratios between friction and inertial 
forces, and between the strengths of the harmonic and 
on-site potentials, respectively. F = Wa/v m . Atkinson 



and Cabrera 11 ll conjectured that only two branches of 
kinks are stable for 

• a branch of static kinks for values of the control 
parameter \F\ below an static threshold F CS (D), 

• a branch of traveling kinks for |f | above a dynamic 
threshold F c d(j, D) < F CS (D), with speeds c larger 
than a minimum speed c c d("f,D). This family has 



a distinctive feature compared to eventual slower 

n 

waves [10| : the leading edge of the kink is monotone 
whereas the trailing edge may develop oscillations. 

The values F cs and F c d correspond to the static and dj 
namic Peierls stresses of the literature on dislocations 
In the overdamped limit 7 — * 00, F cs = F cc i and stable 
wave fronts can be found with arbitrarily small speeds 

Q. 

In a previous paper |10j, we checked numerically the 
validity of Atkinson and Cabrera's conjecture. This is 
a delicate affair and further analytical work is clearly 
desirable. In fact, most numerical studies of kink prop- 
agation truncate the infinite chain to a finite chain, fix 
some boundary conditions and then use a Runge-Kutta 
solver (or variants) to investigate the dynamics of kink- 
like initial configurations. For instance, Peyrard and 



ll 
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Kruskal 1141 applied this procedure to study kinks in 
the conservative Frenkel-Kontorova model, including fric- 
tion near the ends of the truncated chain in an attempt 
to avoid reflections. On the other hand, our analytical 

n 

work |UJ shows that traveling kinks oscillate with almost 
uniform amplitude even for small damping. Then, arti- 
ficial boundary conditions and time discretization may 
greatly distort their shape and dynamics. In fact, using 
Runge-Kutta methods to solve equation J3J with con- 
stant boundary conditions generates reflections at the 
boundary, as shown in Figure ^ after a waiting time 
depending on the size of the lattice. Such oscillations 
end up distorting the right tail and may completely alter 



the shape of the kink giving rise to a complex oscillatory 
pattern. 

A good way to avoid the spurious effects of inappropi- 
ate boundary conditions is to recast as an integral 
equation. Integral reformulations provide an analytical 
expression for the solutions of (0 which we use to de- 
velop new numerical algorithms. Spurious pinning and 
spurious oscillations are suppressed. The introduction 
of these numerical methods based on integral reformula- 
tions of @ is one of our main results. 

The main analytical results of this paper concern the 
nonlinear stability of stationary and traveling wave fronts 
in chains of oscillators. Besides leading to good numer- 
ical methods, we have also used the integral equation 
formulation to investigate the nonlinear stability of wave 
front patterns. We provide a criterion to decide whether 
certain kink-like initial configurations evolve into stable 
wave front patterns. In discrete overdamped models the 
nonlinear stability of traveling wave fronts follows from 
comparison principles. This strategy was applied to the 
study of domain walls in discrete drift-diffusion models 

n 

for semiconductor supcrlattices in |15| . 

Common belief is that comparison principles do not 
hold in models with inertia. This belief is wrong. How 
can we asses the stability of traveling wave fronts in such 
models? For large damping, we can directly compare 
solutions of using its equivalent formulation as an 
integral equation thanks to the positivity of the Green 
functions. As the damping decreases, we can ignore the 



3 

oscillatory tails of the fronts and compare the monotone 
leading edges of the solutions, which drive their motion. 
The process of comparing solutions is technically more 
complex than in the overdamped case because the Green 
functions change sign, and the fronts have oscillatory 
wakes. Summarizing, there are two key ingredients for 
stability. First, the leading edges of the fronts have to 
be monotone. Second, the Green functions of the linear 
problem must be positive for an initial time interval, of 
duration comparable to the time the front needs to ad- 
vance one position. This restricts the possible values of 
the propagation speed for small damping: only fast kinks 
are shown to be stable. Our methods are quite general 
and can be extended to Frenkel-Kontorova models with 
smooth sources at the cost of technical complications. 



The paper is organized as follows. In section 2 we in- 
troduce a numerical algorithm and discuss the stability 
of static kinks. The stability theory for traveling kinks 
is presented in Section 3. In section 4 we discuss the 
role of oscillating Green functions in the appearance of 
static and dynamic thresholds due to coexistence of sta- 
ble static and traveling waves. In section 5 we briefly 
comment on extensions to oscillator chains with smooth 
cubic sources. Section 6 contains our conclusions. Basic 
details on the pertinent Green functions are recalled in 
Appendices A and B. Proofs of our main stability results 
can be found in Appendices C and D. 
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FIG. 1: Trajectory u n (t) computed by solving a truncated 
system of differential equations (dashed) and by integral ex- 
pressions (solid) for 7 = 0.02, D = 4, F = 0.1, n = -70. 



FIG. 2: Convergence to a static kink s n when D = 4, 7 = 10 
and F = 0.05: (a) asterisks u n (t), (b) asterisks 

u' n (t), squares s' n = 0. 



II. STATIC KINKS 

The stationary wave fronts s n for J5J increase from 
s_oo = F to Sqo = 1 + F and solve the second order 
difference equation: 

D(s n+1 - 2s n + s n _ x ) - s n + H(s n - ~) + F = (4) 

in which is the Heaviside unit step function. These 

fronts are translation invariant. We fix their position by 
setting so < \ < s\. Then, s n — F + ar n for n < and 
s n = 1 + F — br~ n for n > 1, where r = 2D+1 +^ 4D+1 -. 
Inserting these formulas in equation for n — and 
n = 1, we find a and 6. Our construction of the stationary 
fronts s n is consistent with the restriction sq < ^ < si 
when |F| < F CS (D). Figure [3 (a) shows a static wave 
front for D = 4 and i 7, = 0.05. As D grows, the number 
of points in the transition layer between the constants 
increases. 



A. Stability 

A stationary wave front s n is stable for the dynamics 
(0) when chains initially close to s n remain near s n for all 
t > 0, as shown in Figure [3 The initial states chosen in 
this figure are u„ = F + <5„ when n < 0, = 1 + F + 5 n 
when n > 1 and = 5*. Both and (5° are small 
random perturbations. 

To find the stable profiles, we proceed as follows. Let 
and u\ be the initial position and velocity of the chain. 
In terms of Green functions calculated in Appendix A, 
u n {t) is given by l|AT2)l with f k (t) = F + H(u k (t) - ±) : 

«»(*) = E*^*(*K + E fc Gk(tK 

+ f Y, k G nk {t-z)H(u k {z)-\)dz (5) 
+FJoE k G° nk (t-z)dz. 



If initially u? < ± for k < 0, ul > § for k > 1: 



«n(*)=£ fc G° fe (t)4 + Gi fc (t)u° + 



(6) 



Jo E fe>0 G° nfc (^ - *)d* + F /J G° nfc (t - z)dz 
as long as tt&(t) < | when fc < 0, tifc(i) > ^ when fc > 1. 
For \F\ < F CS (D), the static wave front s n with sq < ^ < 
si is a solution of equation (0 that satisfies: 



+F!oJ2 k G° nk (t-z)dz 
for all < > 0. Subtracting J7J from |JB}, we obtain: 



(7) 



«n(i) 



This expression holds for t > provided w„(t) — ^ does 
not change sign for any n and t > 0. For which profiles 
u n (t) is this true? Let us select the initial state of the 
chain in the set: 



(9) 



M< J RMin(l,^) 

with R — Min(|s — ||, Si — |) and Co, Ci to be defined 
below. For 7 > 0, we show in Appendix B that \G^ k \ < 
C e-^, \G nk \ < C x e' r > 1 with 7 > 0. This boundedness 
property of the Green functions and © yield: 



K(t)-s n \ < (Co + d)e-*M 



(10) 



Then, \u n (t) — s n \ < R and u n (t) — h cannot change sign 
for any t > 0. Moreover, u n (t) — > s„ as f — > 0. 

In summary, the static kinks are exponentially and 
asymptotically stable in the damped case. Their basin 
of attraction includes all initial configurations and 
selected according to In the conservative case, the 



static kinks are merely stable, but not asymptotically sta- 
ble, because the previous argument with 7 = 0, rj = 0, 
Co = C\ = 1 only yields \u n (t) — s n \ < 2M for all times. 

In the continuum limit D — * 00, the number of points 
in the transition layer between constants increases and 
the distance between points decreases. Then, so and s\ 
tend to \ and the set of states © attracted by s n shrinks 
as D grows. It becomes more likely that initial kinks in 
the chain propagate for a while and finally become pinned 
at some shifted static kink v n = s n+ i, v_\ <\< 

B. Numerical algorithm 

Formula JSJ can be used to compute numerically the 
dynamics of the chain. However, the computational cost 
is high, due to the integral terms and the Green functions. 
In this section, we exploit the static front solutions s n to 
reduce the cost and derive new formulae for u n (t) which 
clarify the dynamics of the chain. 

We will focus on initial kink-like initial states that gen- 
erate ordered dynamics: u k (t) — I changes sign in an 
ordered way as the kink advances. Once the kink has 
passed, the configuration of the chain is close to a shifted 
static kink. That is why we use static kinks to obtain sim- 
plified expressions for u n (t). For instance, let us choose 
a piecewise constant initial profile w° = F for n < and 
u n = 1 + F for n > 1, with u n = 0. For F > 0, Figure El 
shows that — h change sign at time t k , k = 0, 1, 

with to < t% < ... < t k < .... Eventually, the kink may 
get pinned at some static configuration and this process 




FIG. 3: Trajectories u n (t), n 
7 = 0.4, F = 0.15. 

stops at some k. We then use a slightly modified version 
of the integral expression Q for the static wave fronts to 
successively remove the integral terms in JSJ and obtain 
simple formulae for u n (t) similar to ©. In this way, we 
find a relatively cheap algorithm for the computation of 
U n (t). 

Let us describe the algorithm for F > and an initial 
step- like state with Uq < | < ii®, as in Figure |3] 
We must distinguish two cases: < F < F CS (D) and 
F > F CS (D). 



1. Case < F < F CB {D) 

In this case, the stationary wave fronts can be used 
to generate a faster algorithm for obtaining u n (t). we 
remove the integrals in by using the static wave front 
solution of J2J), s n , such that sq < i < s±. 



(12) 



Initial stage. Formula J3J allows to compute u n (t) up 
to the time to at which uo(t) — | changes sign. For t > to 
we compute u n (t) using as initial data u n (to) and u' n (to) 
at to. The later is obtained differentiating JHJ: 

<(*o)=£4^^+^(<-*4 (id 

For t < t < ti, equation (JSJ becomes: 

«n (*) = Ek k (t - to)u' k (t ) + G\ k (t - to)u k (to) 

+ It £*>_! G° fe (t - z)dz+Ffl £ fc G° fe (i - z)d 2 
Now, u_i(to) < \ < uo(to) and we must use the shifted 
stationary solution v n = s n +i, that satisfies u_i < 5 < 
t>o- Observing that s„+i solves (|2J) with initial data 
s n +i, at time io we obtain the formula: 

sn+i = Efe G r \ fc (t - t )**+i + F /* £ fe G« fc (i - z)<fc 

+ / t t E fc >-iG° fc (i-z)^ 
Subtracting (|13|l from (|12|l we find: 

u»(f) = s„+i + E fc G° fc (i - *o)«' fc (*o) 

+ E*Gi fc (t-to)(u fc (to)-a fc+ i) ) 
up to the time t\ at which u-i(t) — i changes sign. 

Generic step. Once we have computed the time ti at 
which ui (t) — 5 changes sign, we calculate the new initial 
data u n {ti) and 

u»(*0 = + E fe G£ fc fa - 

+ E fc G^(ti-*i-i)(«fc(*i)-*fc+i). 



(13) 



(14) 



(15) 



<(*«) = E*^(*i-*i-iK(*0+ 



(16) 



Then the evolution of the chain for t > ti is given by the 
formula: 

Un(t) = S n+l+1 + J2 k G° nk (t - t,)u' k (tl) + 



(17) 



Efc G nfc0 - U)iu k (ti) - Sfc+J+l), 



until either u_(; +1 )(t) — \ or U-i(t) — | change sign. If 
— 5 changes its sign at a time we start a 
new step using to compute u n (t). If U-i(t) — I 

reverses its sign at a time tj+i, we start a new step using 
Sn+2-i to compute u n (t). 



2. CaseF>F ca {D) 

In this case, it is convenient to remove the integral in 
Eq. (JSJ) by using as s n the static wave front solution of (J2J) 
corresponding to an applied force F = F CS (D), and such 
that So < \ < s\. Recall that there are no stationary 
wave fronts for F > F cs . 

Initial stage. Subtracting (JJJ at F CS (D) from JSJ we 
find: 



unit) =s n + £ fe G nk {t)u\ + £ fc Gl k {t){v? k - s k ) 
+(F-F cs )J^ k G° nk (t-z)dz 



The remaining integral term can be removed by observing 
that 1 is a solution of J2J with F = and initial data 
u n {0) = 1, <(0) = 0: 

1 = E G nk (*) + / E G ^( f - ( 19 ) 

h Ja k 

Multiplying (|19|l by (F — F cs ) and inserting the result in 
(|18|l we obtain: 



^< + ^(ul-s k -F+F cs j 



up to the time to at which uo(t) — \ changes sign. For 
t > t : 

u n (t) = s n+1 + (F — F cs ) + J2 k G° nk (t - t )u' k (t ) 

+ E k G] lk {t - to)(« fc (to) - s k+1 — (F — F cs )), (21) 

<(^o)=Efe 

up to the time t\ at which u-\(t) — i changes sign. 

Generic step. Similar to the generic step for F < F cs 
but replacing s n by s n + (F ~ F cs ). 

3. Numerical implementation 

We will use I|17f) - (I16|I and l|21|) to study the dynamics of 
the chain in Section IlIII Due to translational invariancc 
Glk = G l-kfi and G lk = g «-m- To calculate u n (t), we 
only need to compute G" (t), G} lQ (t) for a time interval 
[0,T], T < max/|t i+ i - U\ and for \n\ < N, where N 
is sufficiently large. We calculate the integrals G° (t), 



(18) GUt) 



dt 



dt 



by means of the Simpson rule, 



choosing a step smaller than the period of the oscillatory 
factors. The value N is selected so as to make the error 
introduced by the truncated series J2\n-k\<N sufficiently 
small. This is possible because the Green functions and 
their derivatives decay as \n — k\ grows. 

A more general version of our algorithm will be pre- 
sented elsewhere 



i ve i 
If]. 



III. STABILITY OF TRAVELING KINKS 



Un(t) = S n + (F~F CS )+ J^k G° nk (t)ul 

+ E k G 1 nk (t)(u k -s k -(F-F cs )) 



(20) 



In this section we introduce a strategy to study the 
stability of traveling wave fronts in J3J). 
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Traveling wave fronts are constructed by inserting 
w n (t) = w(n — ct) in (J2J to produce a nonlinear eigen- 
value problem for the profile w(x) and the speed c. As- 
suming w(x) < i for x < and w(x) > h for x > 0, 
the problem becomes linear. The wave profiles are com- 
puted as contour integrals, imposing w(0) — \ to find 

nn 

a relationship between c and F |10t The law F(c) 
and the shape of the wave profiles are controlled by the 
poles contributing to the contour integrals. The rele- 
vant poles depend on the strength of the damping. For 
large damping, we have complex poles with large imag- 
inary parts. The dependence law F(c) is monotonically 
increasing and the wave profiles are monotone. For small 
damping, poles with small imaginary parts become rele- 
vant, in increasing number as the speed c decreases. The 
function F(c) oscillates for small speeds. Different os- 
cillatory wave profiles with different speeds may coexist 
for the same F . At zero damping, those poles become 
real and the wave profiles develop non decaying oscilla- 
tions. For some ranges of speeds, the waves constructed 
in this way violate the restriction w(x) < ^ for x < and 
w(x) > \ for x > 0. Those ranges should be investigated 
with a modified technique allowing for a finite number of 
turning points. 



Complex variable methods yield families of explicit 
wave solutions but give no information on their stabil- 
ity. Numerical tests and physical context Q suggest 
the stability of traveling kinks that have monotone lead- 
ing edges and large enough speeds. Figures depict 



the wave profiles for decreasing 7. We now confirm that 
these wave fronts are stable. The travelling wave w n (t) is 
stable for the dynamics of the chain when the solutions 
Unit) of J5J) remain near w n (t) for all t > if the initial 
states u\ are chosen near w„(0), w' n (0). Controlling 
the evolution of u n (t) is more or less difficult depending 
on the properties of the Green functions. We distinguish 
two cases: positive Green functions (large damping) and 
oscillatory Green functions (small damping). 

A. Strong damping 

For large damping 7 2 3> 4, we know that the wave 
front profiles are monotonically increasing and that the 
Green functions are positive and decay exponentially in 
time (cf. Appendix B). The main result of this section is 
the following stability theorem, whose proof can be found 
in Appendix C: 

Theorem. Let us select the wave front profile so that 
w n (t) = w(n — ct — i), with c < 0, and F > 0. // we 
choose the initial states for u\ and u^, in the set: 

1 



w n {-T) < < < w„(r), 0<t< 



2\c\ 



ther 



WJ-t) - <| < u° - w„(-t) 

\ W 'n( T ) ~ U i\ < w n(r) ~ < 



U n (t - T ) < u n(t) < W n (t + t) 



(22) 



(23) 



(24) 



for all n and t > 0. 

In other words, if the initial oscillator configuration is 
sandwiched between two wave front profiles with differ- 
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ent phase shifts, w n (— t) and w n {r), (with a sufficiently 
small r), then the oscillator chain remains trapped be- 
tween the two shifted profiles w n (t — r) and w n (t — r) 
forever, provided {u^ — w' n (0)\ is sufficiently small. This 
implies the dynamical stability of the wave. The more 
involved argument explained in Section [ill Bl for conser- 
vative dynamics can be used to prove that the wave fronts 
are also asymptotically stable. 

Furthermore, the basin of attraction of a particular 
traveling wave is larger than 11' I'll - (|23|l . as shown in 
Figures BISI for F > F CS (D). The initial oscillator con- 
figuration in this figure is a step function, u® n = F for 
n < and it„ = 1 + F for n > 1, with a superimposed 
small random disturbance. The initial velocity profile 
fluctuates randomly about zero with a small amplitude. 
After an initial transient, the trajectories get trapped be- 
tween advanced wave fronts, w n (t + T), and delayed wave 
fronts, w n (t — t). Moreover, they converge to a shifted 
wave front, w n (t + a), as t — > oo. 

B. Conservative dynamics 

For small damping 7 2 -C 4, we know that the kink 
profiles develop oscillations in the trailing edge (see Fig- 
ure [JjJ and that the Green functions oscillate and change 
sign (cf. Appendix B). However, G° fe (f) and G\ k (t) are 
positive for < t < T* = T(j, D). This critical time T* 
plays a key role for the stable propagation of waves. We 
will show in this section that kinks are stable provided 
|c| > 7^-. Our argument does not say anything about the 



1.2 



0.8 



0.4 




FIG. 4: For D = 4, 7 = 2.2, and F = 0.2: (a) Compared time 
evolution of wj n (t+r) (dot-dashed line) w n (t — r) (dashed line) 
and u„(t) (solid line) when n = 0, —1, —2, (b) Compared 
profiles u n (T) (circles), w„(T ± t) (asterisks), (c) Compared 
time evolution of w' n (t±T) (dot-dashed and dashed lines) and 
u' n (t) (solid line). 
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FIG. 5: Same as Figure H when F = 0.45. 

stability of kinks with lesser speeds. Moreover, T* — > 
and our lower bound on the wave front velocity tends to 
infinity, in the continuum limit. 

We show in Appendix B that a rough estimate for T* 
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is provided by 



For 7 = and D 



% /4(l+4_D)- 7 2 

as chosen in our Figures 10-0, T* > 1. Then, kinks 



Q, 



with |c| > 1 are stable. In [llj and [lfj, stability was 
conjectured for speeds larger than the last minimum of 
F(c), which is attained at c c d ~ 0.74. 

For small or zero damping we cannot use the previ- 
ous comparison arguments because the trailing edge of 
the traveling wave front oscillates and monotonicity does 
not hold there. If we look at the traveling wave front 
profiles, it becomes clear that we should compare the 
monotone leading edges of the fronts. Figure [7| (a)-(b) 
depicts the trajectories w n (t) and their time derivatives 
w' n (t) for a particular traveling wave front. We observe 
that w n -x(t) < w n (t) < w n+ i(t) and w' n _ 1 (t) < w' n (t) < 
w' n+1 (t) up to a certain time. Figure [7| (c) shows the 
initial configurations for w n (0) and the shifted waves 
w n (— t), w n (r). uVj(Q) is sandwiched between w n (— r) 
and w n (r) up to a point no- Figure 0(d) depicts the ini- 
tial velocity profiles w' n (0), w' n {—T) and «4(t). w' n (0) is 
sandwiched between w' n {— t) and w' n {r) up to a point n%. 
no and n\ mark the onset of the oscillatory tails. In gen- 
eral, < no < rii. As the wave advances, the ranges of n 
for which w n {t — t) < w n (t) = w(n — ct — |) < w n (t + r) 
change with t. 

The main result of this section is the following stability 
theorem, whose proof can be found in Appendix D: 

Theorem. Let us select the wave front profile so that 
w n{t) — iv (n — ct — i), with c < and F > 0. Let T* 
be the maximum time up to which the Green functions 



G^ k (t) and G^ k (t) remain positive. We assume that the 
speed \c\ > and choose the initial states for u x n 
and u„, in the set: 



vj„(-t) < u° < W„(t), n<n < t < ^j- 



(25) 



w' n {-T) <u\< w' n (r), n<m 

E„K-»n(0)|< £ , £„|<-<(0)|<e. (26) 

for e > small and r < e. Then, we can find an increas- 
ing sequence of times tk, k = Q, 1, with t_i = 0, such 
that: 



Wn(t — t) < u n (t) < w n (t + t) , n < no — k 



(27) 



w' n (t- t) < u' n (t) < w' n (t + T), n<n x -k 

for tk-i <t<tk- Furthermore, for t > and any n, we 
have: 

M*) - w n {t)\ < K k (t)H - Mo)\ 
k 

+ J2\ G no(t)\\ul-w' k (0)\ + C(t) (28) 

k 

C(t)= Y, / G° nk (t-z)dz, 

Ai<0,t>T fc -r Tfc ~ T 

in which T k — j^j- + ■ Thus, the traveling wave front is 
stable when 7 = or asymptotically stable when 7 > 0. 



Let us clarify the meaning of (|28p. For 7 > 
0, the sums £fel G nJfc(*)l> Efcl G °fcWI deca y exponen- 
tially with time. For small t, the function \C(t)\ ~ 

2r Sfe<o,t>T fc -r l^nfc(* _7 fe)l- Tnis sum is nnite an d de- 
cays with time. This explains our asymptotical stability 
claim. When 7 = 0, the sums £ fe | G?^ fc (*) 1 1 rtg - w k (0)\, 
Sfc l^nfc(*)ll u fc — w k{0)\ are bounded by a constant times 
Uax k \u o k -Wk(0)\+Max k \u° k -w k (0)\. The function \C(t)\ 
is bounded by a constant times r and is made small be 
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FIG. 6: Trajectories u n (t) when L> = 4: (a) 7 = 0.2, F = 0.2, 
(b) 7 = 0.2, F = 0.45, (c) 7 = 0.1, F = 0.2, (d) 7 = 0.1, 
F = 0.45. 

choosing r small. This explains our stability claim in the 
conservative case. 

The inequalities (|27|l tell us that the leading edge of 
the propagating kink is sandwiched between the leading 
edges of the shifted traveling wave fronts w n (t + t) and 
w n {t — t). As the kink u n (t) advances, the times tk at 
which U-kif) — h changes sign are bounded by the times 
at which the advanced and delayed wave fronts cross i: 
Tfe — t < tk < Tk + t. This fact is the key for obtaining 
the stability bound . 

IV. COEXISTENCE 



The results in Section fill Bl indicate that stable static 
and traveling kinks may coexist. The only restriction 
on the traveling kinks is the monotonicity of the leading 
edge and a low bound on the speed. These conditions 



FIG. 7: (a) Trajectories w n -i(i) (dash-dotted), w n (i) (solid), 
w„+i(t) (dashed); (b) same for w'„_ 1 (t), w' n (t), w n +i(t); (c) 
initial configurations for w n (r) (circles), w n (0) (asterisks), 
w n (— t) (squares), the vertical line defines no; (d) same for 
w' n {r), u4(0), w' n (— r), the vertical lines defines m. 

are satisfied by traveling wave fronts for a range of forces 
in which static wave fronts also exist. We show in this 
section how oscillating Green functions may force initial 
kink-like configurations (which would be pinned for large 
damping) to evolve into a traveling wave fronts provided 
the damping is small enough. 

We fix F < F cs and select the static kink s n con- 
structed in Section [Q] for J2J with so < | < Si. Let 
the initial condition for @ be a piecewise constant pro- 
file: u° = F for n < Q, u° n = 1 + F for n > 1 and u\ = 0. 
Let T* the maximum time up to which G° nk and G\ k 
remain positive. 

As long as u n (t) — i does not change sign for any n, 
u n (t) is given by formula Q in Section ^1 We have 
E fc Gl k {t)u\ > for t < T*. Initially, G] lk (t) is concen- 



trated at k = n and the sign of J2k Gnk(t)( u k ~ s fc) is 
decided by the sign of — s n . If w° > s„, u n (t) > s n - 
In our case, this is true for n > 1. If u„ < s n , u n (t) 
increases towards s„ as ^2 k G\ k (t){u° k — Sk) decays. By 
our choice of the initial state, uo(t) grows faster than the 
other components u n (t), n < 0. 

Now there are two possibilities depending on the value 
of the damping coefficients. For large damping, G\ k (t) 
is positive for all times and G° nk decays fast. Then Uo(t) 
cannot surpass sq. This initial data is pinned for large 
damping. 

For small values of the damping, G\ k (t) changes sign. 
Then uo(t) given by Q) may surpass sq and i since the 
term J2k ^nk(^)( u k — s k) becomes positive for t > T* . 
This process can be iterated to get a stably propagating 
wave, see Figure ^ A prediction for the speed is found 
in this way: it is the reciprocal of the time that sq + 
Gl (t)(F — s ) needs to reach i. 

V. MORE GENERAL POTENTIALS 

We have focused our study on periodic piecewise 
parabolic potentials V(u) — \u\ < i. For these po- 
tentials, families of static and traveling wave fronts can be 
constructed analytically. Schmidt 17j and later authors 
|lJ. found exact monotone wave fronts of conserva- 
tive systems by constructing models with nonlinearities 
such that the desired wave fronts were solutions of the 
models. For damped Frenkcl-Kontorova or quartic dou- 
ble well potentials, stably propagating wave fronts have 
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been found numerically 10]. Numerical studies of kink 
propagation in the conservative Frenkel-Kontorova model 
were carried out in 

The stability of propagating kinks in these models can 
be studied adapting the methods developed in this paper, 
but the analysis is technically more complicated [lj| • For 



:onsc 

Q. 



instance, taking V{u) — %- for |u| < |, 
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2 ^ 2 , 4 2 

\u — 1| < 2 w e get a continuous piecewise linear source 



for 



g(s) 



s s < 1/2 



I i<«<3/2 



s-1 s>3/2 
The arguments in Sections llII Al and lHI Bl can be adapted 
by including new terms Jo£| G ° nk (t - z)(2u k (z) - 
\)dz in the integral expressions (|C1I) - (|C2|) and using that 
the function h(s) = 2s — 1/2 is increasing in 1/2 < 
s < 3/2. Similarly, for a Frenkel-Kontorova potential, 
we write g(s) = —as + (sin(s) + as), a > 0. Then, we 
find the integral expression 1A12|> with a nonlinear source 
fk = — sin(ufc) + auk + F, using the Green functions for 
the linear operator u'^+ju' n — D(u n+ i — 2u n +u n ^i)+au n . 
The parameter a is chosen to ensure adequate monotonic- 



equa 



ity properties for h(s) = — sin(s) + as [16|. 

VI. CONCLUSIONS 

We have developed a nonlinear stability theory 
for wave fronts in conservative and damped Frenkel- 
Kontorova models with piecewise linear sources based on 
integral formulations. Our results provide an analyti- 
cal basis for the distinction between static and dynamic 
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Peierls stresses, which arise as thresholds for the exis- 
tence of stable static and traveling wave fronts. With 
little or zero damping, stable propagation of fronts is pos- 
sible when their speeds surpass a critical value. The cor- 
responding wave front profiles have a monotone leading 
edge, and, possibly, an oscillatory wake. Wave fronts can 
be oscillatory and yet stable. Whether slow wave fronts 
showing oscillations in the leading and trailing edges are 
stable remains an open question |ic| . It is remarkable 
that high order quasicontinuum approximations such as 

n n 

those by Rosenau 20] or by Boussinesq 21] have wave 
solutions comparable to the fast waves of the discrete 
conservative model [l2| . 

Together with the stability theory we have presented 
an algorithm for the numerical computation of the dy- 
namics of kinks. Our scheme has good stability proper- 
ties and avoids distortions originated by artificial bound- 
ary conditions and time discretization. 
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APPENDIX A: GREEN FUNCTIONS 

We want to find an integral representation of the solu- 
tion of the problem: 

u'n + ju' n = D(u n+1 - 2u n + u n -\) -u n + f n (Al) 
Un(0)=«°, <(0) = u* (A2) 

with 7 > 0, D > 0. Firstly, we get rid of the difference 
operator by using the generating functions p(9, t), f(9, t): 

p(9,t) = J2 u n(t)e- m6 , Mt) =Y,f-{t)e- ine (A3) 

n n 

Differentiating p with respect to t and using I|A1J) . we see 
that p solves the ordinary differential equation: 

p"(6, t) + 7 p'(0, t) + u;(9) 2 p(9, t) = f(9, t) (A4) 

where oj(9) 2 = 1 + AD sin 2 (|) and the obvious initial 
conditions for p(8,t) follow from those for u n (t). 

The solution p depends on the roots of the polynomial 
r 2 + 7 r + uj(6f = 0. When £ > u 2 {6), 

e M«)* r+ (0)-e r +W t r_(0) 



p(M)=p(0,O). 



+p'(o, o) 



r+(0)-r_(0) 

e r+(0)t _ e r-(e)t 



r + {6)-r_{e) 



* e r+(e)(t-s) _ e r_(fl)(t-«) 

o r+(fl)-r_(0) 



-f(9,s)ds (A5) 



with r±{9) 



-1±Ja{6) 



< and a(6) = 7 2 - Alo{9) 2 



When ^j- < oj 2 {9), the roots are complex: 



p(M)=p(MK 



cos 



. , s , 7 sin (J(0)t) 



+p'(0,0) e~ 



^ t sin (I{9)t) 



e 2 



sin (1(9) (t-sj) 
1(9) 



1(9) 
f(6,s)ds (A6) 
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where 1(9) = ^ ° W . When £ = w 2 ((9), 



p(9, t) = p(6, 0)e^*(l + |t) + 0) e^*i 



+ / (t-s)e^(^ s) /(^,s)rfs (A7) 



The solution u n (t) of ijAlf) is recovered from the defi- 
nition l[A3j) : 



/ £e m6 p(9,t). 



Here, p(6>, t) is defined by (|A5|) when 9 £ I± 



h = {9e[-7:,n}\^>u J 2 (9)}, 



(A8) 



(A9) 



by HA6|> when 9 e h 



h = {9e [-7T,7r] I ^ <w 2 (6»)}. (A10) 



and by l|A7[) when e J 3 



P = {0G[-7T,7r]|^=a; 2 (0)}. (All) 



Notice that h = P = if 7 2 < 4 and 7 2 = P = if 
7 2 > 4(1 + AD). P ^ only if 4(1 + AD) < 7 2 < 4 and 
it then consists of two points. 

Formula l|A8(l can be rewritten as: 

u n {t) = £ [G£*(*K(0) + G* fc (t)u fc (0)] 
fc 

+ fY. G nu{t-s)h{s)ds (A12) 
Jo fc 



where 



G° fc (*) = r w £e^)V(M) 
G^(*)=r w Se <(B -*V(*,t) 



(A13) 



with 



(«)*_„'•- (»)* 



r+(0)-r_(0) ' 

e^H, 9eP 
— 7(e) — 



G Ji 



(A14) 



G/ 2 



/(M)= < 



6 r -' 9 "r + (g)-e'-+" l 'V-(9) 
r + (S)-r_(9) ' 

e^*(l + ^), 6» eP 



G/i 



(A15) 



cos (/(0)t) 



7 sin (l(e)tl 



For conservative chains, 7 = 0, G^. 



dt 



G/ 2 



and: 



r/0 P i(n-k)6 

Gnk(t) = J n h^w s[n{u}m ' (A16) 

Green functions for hamiltonian chains were studied in 
and earlier in . For overdamped chains, they were 



computed in 
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APPENDIX B: PROPERTIES OF THE GREEN 
FUNCTIONS 

The Green functions for (|A1|) - (|A2|) have three relevant 
properties: they decay in time, they decay as \n— k\ — * oo 
and are positive for some time. The property of spatial 
decay follows from integration by parts in 1A13|) : 
r° (t\ — f n de r i(n-k)e d l g°(e,t) 

^nfcW — ii(n-k) 1 J-7T 27r e W N 

(B1J 

r 1 (+\ — r d6 „i(n-k)e d'g 1 (8,t) 

Xj vk\ b ) — i'(Ti-fc) 1 J-tt 27i- e d? 7 

when ?i 7^ fc. An immediate consequence is that 
E fe \G°nk(t)\ P ^d £ fe |G^MI P are finite for any p > 1. 
Therefore, we may obtain decay results as |n| — > co for 
the solutions u n (t) of (|A1 |) -(|A2 |I given by (|A12|I decay 
when the data m° , it* , / n (t) decay. Figures 18131 illustrate 
the spatial decay of G° n (t) = G° nQ (t) and G l n (t) = G£ (t). 
Notice that, initially, both are concentrated about n = 0. 

Time decay and positivity depend on the strength of 
the damping. Let us start by the strongly damped case: 
7 2 > 4(1+4D). The Green functions are given by l|A13|l - 
(lA15fl with I 2 = P = 0: 
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^nfc(^) an d ^nfe(^) are positive- The roots r±(8) 
being even with respect to 9, both G° lk (t) and 
G] lk (t) are real and e l ("- fc ) 6 ' can be replaced by 
cos((n - fc)0). The kernels g (9,t) = e ^ +(e) Z e r ^ { g) 
and.gi(6»,t) = ^gp^gy r ~ W are even, pos- 

itive, reach their maximum values at = and de- 
cay as 9 increases to 7r. The dominant contribution 
to the integrals l|A13() comes from a neighbourhood 
centered at 9 — 0, where the oscillatory multiplier 
cos((n — k)9) is positive. Thus, both G^ k (t) and 
G] lk {t) are positive. Figure 0D illustrates their evo- 
lution as time grows. Notice the resemblance with 
the time evolution of heat kernels. 

(*nk W an d G* fc W are bounded uniformly in n, k 
by decaying exponentials in time: 

r + (0)t_ r_(ir)i 



|G^)| < 



(B2) 



v / 7 2 -4(l+4_D) 

We come now to intermediate damping 4 < j 2 < 
4(1 + AD). In this case both I\ and I2 are non empty. 
The piecewise defined kernels go and 171 are still even, 
take the largest values near zero (in Ji) and the small- 
est near ir (in 12). The dominant contribution to G a nk (t) 
and Git,(t) comes thus from ii and is positive. This is 
helped by the fact that the contribution coming from I2 is 
initially positive and the factor e - ^* in the oscillatory re- 
gion I2 decays faster than the factor e r +^ t in the positive 
region I\. Therefore, G° fc (i) and G^ fe (£) are essentially 
positive in this intermediate regime, see Figure |5J This 
means that their large components are positive, despite 
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FIG. 8: Time evolution of the Green functions for D = 4 and 
7 = 10: (a) t = 0, (b) t = 0.5, (c) t = 5; (d) t = 0, (e) t = 0.5, 
(f) t = 5. 
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FIG. 9: Time evolution of the Green functions for D = 4 and 
7 = 2.2: (a) and (d) t = 0.5, (b) and (e) t = 2.5, (c) and (f) 
t = 5. 

the appearance of some negligible negative components. 
They can be roughly bounded by: 

\G° nk (t)\ < C e r+(0)t , \Gl k {t)\ < G ie ^°>' (B3) 

We address finally the weakly damped problems with 
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FIG. 10: Time evolution of the Green functions for D = 4 
and 7=1: (a) t = 0, (b) t = 0.5, (c) t = 2.5; (d) t = 0.5, (e) 
t = 2.5, (f) t = 5. 
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FIG. 11: Time evolution of the Green functions for D = 4 
and 7 = 0: (a) and (d) t = 0.05, (b) and (e) t = 0.5, (c) and 
(f) t = 2.5. 

7 2 < 4. In this case, I 1 =P = H). G" nk {t) and G* fc (t) are 
no longer globally positive. However, the kernels go(6,t) 
and gi(0,t) are positive for \t\ < —, = = 2T 



and lil < 



T, respectively. That means 



1 /4(l+4D)-7 2 

that G a nk (t) > and G^t) > for t in those intervals. 



They remain essentially positive in a larger interval. The 
kernels go(0,t) and g±(0,t) become negative for 6 near 
7T and remain positive for small 9. This is enough for 
the relevant values of G x nk (t) to remain positive up to a 
critical time T* , often larger than T. We can get uniform 
bounds in time: 



\G° nk (t)\ < 



2 — ~-Lt 

- — e 2* 



(B4) 



\Gl lk (t)\<(l + 7 ^)e-¥ 



The same positivity properties and bounds are shared 
by the Green functions in the conservative case 7 = 0. 
Figures EH an( l [HI illustrate the time evolution of the 
Green functions. A detailed study of the decay properties 
with respect to n and t for conservative problems can be 
found in 



APPENDIX C: STABILITY OF TRAVELING 
WAVE FRONTS FOR STRONG DAMPING 



We now prove the stability theorem of Section UlI Al for 
strong damping. The key idea of the proof is suggested 
by formula JSJ . When we solve J2J) starting from different 
step-like initial states, we observe three types of terms in 
0. The second and the third are increasing functions 
of the step-like configurations. The fourth term does not 
depend on the initial configuration. The first one can 
be made small by choosing a small velocity profile. Our 
proof proceeds in two steps. First, we establish a few 
properties of the traveling wave fronts. Second, we prove 
the stability bound ipijl . 
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Step 1: Basic properties of the traveling waves. For ev- 
ery k, we know that w' k (t) > 0. Thus, each w k {t) crosses 
\ at a definite time t k . Recall that we have selected the 
unique wave profile w(x) satisfying w(0) = ^. Therefore, 
w- k {t) - \ changes sign at time T k = ^ + ^ > 0, 
k = 0,1,2,.... For the shifted waves Wk(t + r) and 
Wk{t — t) the changes of sign take place at the shifted 
times T k = T k — t and T k = T k + r. 

The waves w n (t ± r) solve the integral equation Q 
with initial data w n (±r) and w^(±r). Using the times 
T k , we can rewrite formula (0 in a more explicit form: 

w n (t± T )=J2 [G° fc (t)w' k (±r) + G l nk w k (±t) 

k 

+ F / E G nk(t ~ *)dz + [J2 G ^ ~ Z ^ Z 
Ja h J ° fc>0 

/•Max(t,T±) 

+ E/ ± G° nk (t-z)dz (CI) 

fc<0 T k 

A term J T ± G^ k (t — z)g?2: is added whenever a factor 
iDfc (2 ± r) — I changes sign. 

S^ep 2: Comparing u n (t) and w n (t ± r). During the 
initial stage of the evolution of the chain uo(t) < \ < 
u\{t) and formula © reads: 



Therefore, HClfHC41 imply 



So E fe>0 G° nk (t - z)dz + Fj* J2 k G° nk (t z)dz 
By and the positivity of G* fc (t), 



(C2) 



Tr 1 

k 



pl lk (t)M-r) <YPlk(thl <YG 1 nk {t)Mr)- (C3) 



By 



EkGL(t)[Mr)-ul]+G" nk (t)[w' k (r)-ul] > 



(C4) 



w n (t - t) < u n (t) < w„(t + r), 



(C5) 



for all n and t <T Q + . Recall that T + < T ~ by definition. 
Afterwards, wo(t + t) has crossed \ and f^,+ G° (t — 
> must be added in the expression for w n (t + r). 
The ordering l|H5]) still holds. At time T ~~, w (t - t) 
crosses \. By (|C4|I . «o(i) must cross before, at a time io- 
In this way, we obtain a sequence of times t k at which 
U-fc(t) — |, changes sign satisfying Trj" < t k < T k , fc = 
0,1,2,.... Then, 

Un (t)=£[G° fc (*K + G^(*K] + 
fe 

^ / E G °nk(t ~ z)dz + G °nk(t - z)dz 

J ° fe J ° k>0 

/Max(t,tfc) 
G° nk (t-z)dz (C6) 

and l|C5|l holds for all t. Our stability claim is proved. 

APPENDIX D: STABILITY OF TRAVELING 
WAVE FRONTS FOR CONSERVATIVE 
DYNAMICS 

In this section, we prove the stability theorem of Sec- 
tion 1111 Bl for small or zero damping. The notation is the 
same as in Appendix ICl and the proof is organized in two 
steps. 

Step 1: Initial stage. We compare u n (t) given by (|C2|) 
with the shifted waves w n (t ± r) given by I|C1(I . whereas 
u' n (t) is compared with w' n (t ± r). The time derivatives 
are calculated by differentiating ljCl(l - l|C2|l . Notice that 



Efe G\ k (t) K " M-r)] +GI k (t) [u\ - w' k (-r)] > 



dt 



> for small damping when t < T* . Up to a 
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first critical time Tq , u k (t) — | and w k (t ± r) — ^ keep 
their sign for all k. Therefore, 

w n (t—T)<u n (t)<w n (t+T), n<n , 

(Dl) 

W 'n( t - T )< u 'n( t )<' w 'n( t + T ), n < Til, 

for i < T Q + . Recall that, initially, G\ k and G^ fc , together 
with their derivatives, take on their maximum values for 
k close to n. This fact and l|2"I )l -l|2"5 |l imply: 

w n (t - t) < u n (t) < w n (t + r), n < n - 1, 

(D2) 

<(t - r) < < w' n (t +t), n<m- 1, 

for T+ < t < T+, choosing r < T* — T . This means 
that uo(t) — \ changes sign at a time to such that T + < 
t < Tq - < T x + . We then obtain formula l|C6p for u n (t) 
restricted to t < Ty. By subtracting ijClf) from l|C6|) . we 
find: 

E„ K - w n \(t) < j: n \Gl (t)\ E n K - u>„|(o) 



E„ |G° (*)I £„!<-<! (o) + c(*) 



(D3) 



where 



C(t) = 



— > 



Q G° (t-z)^, t>T+, 



and 



E„ K - <\(t) < En l^(*)l E„ l«n " tWnKO) 



EJ%K*)IEJ<-<l(o) + m 



(D4) 



where 



+ 



0, t < T G , 



for i < T+. Let Si = Max [0lOO ) E„ |^(*)|, S 2 = 
Max[o,oo)Enl^(*)l. ^ = Max [0!Oo) E„l^(i)l and 
5 4 = Max [0)Oo) E„ I^WI- Then, for t < T+, 

E„ M*) - <"n(*)| < (ft + ^ 2 )e + 2r|G° (T )|, 



E„ K(*) " <(*)! < (ft + ^ 4 )e + 2r|^(T )| 



(D5) 



The distances \u n (t) — w n (t)\, \u' n (t) — w' n (t)\ remain of 
order s. In particular, the oscillatory tail of u n (t) for 
n > no is contained in the same band that contains w n (t) 
for n > no. 

Step 2: Generic stage. We iterate Step 1 starting at 
times T ; + , I = 1,2,..., according to the following induc- 
tion procedure. For a fixed T ; + , 1D2|) holds for n < no — I, 
n<m-l, T l t 1 <t < T ; + and: 



\u n (t) - w n (t)\ < (Si + S 2 )e + 2tS, 



(D6) 



E n \<(t)~<(t)\<(S 3 + S 4 )e + 2TS 
holds for t < T+, with S = Max(£ fc |G° fc (T fc )|, 

dG° 

Efe I d ^ fc (^fc) |) - Now, we shall show that these properties 
also hold for 

For T ; + < t < Tj\ v the evolution of w n (t ± r) and 
u n (t) is given by: 

w n (t ± r) = £ fe G° fc (t - T+)«/ fc (T+ ± r)+ 
E,G J \ fc (<-T+K(r+±r)+ 



(D7) 



"n(i)=E & G° fe (t-T+K(T+)+ 

E fc ^ fc (t-T+)u fe (T+)+ 

Ff T +L k G nk {t-z)dz+ 



(D8) 



J T +Efc>-;G£ fe (>-2)cte+ 
Efc<-i Jt+ G nk( f ~ z)H{u k {z) - \)dz. 
Notice that we have taken as initial data the values of 
w n (t ± r) and u n {t) at time T ; + . In this way, formu- 
lae (JD7(1 and i|D8|) only involve the values of the Green 



functions in a short time interval [0, Ar + 2r], Since 
A + 2t > T*, the Green functions are both positive. 
Recall that for this short time interval G x nk and G° k , to- 
gether with their derivatives, take on large values for k 
close to n. We can then use (|D2(1 for n < no — l,n < n\ — l 
at time T+, jD6j and {D7jl-l|D8l) to obtain {D2J for 
?i < n -(7 + 1), ?i < tii-(/ + 1) and T+ < t < T+ r This 
means that U-i(t) — \ changes sign at a time ti such that 
T+ < ti < Tf <T+ +l . We then obtain formula (pi for 
u n (t) restricted to t < T^_ v Subtracting (fCl"]) from (fC6|) 
for t < T+ +1 , we find: 

E„ K - Wn\(t) < (S 1 + S 2 )s + 2r E,< ; K k (T k )\ 
E„ l< - <IW < (S3 + S,)e + 2rj: k<l \^(T k )\ 
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This implies (|D6p for t < T^. We are now ready to 
repeat the process starting at time T^_ v 

Step 3: Conclusion. From Step 2 we obtain a sequence 
of times ti for I = 1,2, with T + < U < , at which 
U-i{t) — i changes sign. In this way, we keep track of the 
times ti at which changes of sign take place and obtain 
formula l|C6() for u n (t) for all t. Subtracting i|Cl(l from 
(IC6I) we find the bound l|28|) on \u n — w n \. 
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